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0^ ■ Abstract 

-«^ ■ We present a new analysis of the ratio e'/e which measures the direct CP vi- 

o ■ olation in K ^^ vrvr decays. We use the 1/Nc expansion within the framework of 

O ; the effective chiral lagrangian for pseudoscalar mesons. The 1/Nc corrections to 






X 



the hadronic matrix elements of all operators are calculated at leading order in the 
chiral expansion. Performing a simple scanning of the input parameters we obtain 
1.5 • 10~^ < e' /e < 31.6 • 10~^. We also investigate, in the chiral limit, the l/N^ 



c 



;h ' corrections to the operator Qq at next-to-leading order in the chiral expansion. We 

find large positive corrections which further enhance e'/e and can bring the standard 
model prediction close to the measured value for central values of the parameters. 
Our result indicates that at the level of the 1/Nc corrections a A/ = 1/2 enhance- 
ment is operative for Qq similar to the one of Qi and Q2 which dominate the CP 
conserving amplitude. 
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1 Introduction 

There are two types of CP violation which appear in the neutral kaon system: direct and 
indirect. Direct CP violation occurs in the amplitudes and will be the subject of this 
paper. Indirect violation occurs in the physical states and is characterized in a phase 
convention independent way by the parameter e. Indirect CP violation has been observed 
and is incorporated in the standard model, as a restriction to the CKM phase. Direct CP 
violation is described by the parameter e' whose predictions require greater attention and 
has been the subject of several investigations. The superweak theory |l| predicts e' je to 
be exactly zero. In the standard model the predictions for the ratio cover a wide range of 
values. Until recently, the experimental evidence for e' /e was inconclusive. While the value 
Re {e' /e) = (23 ±7) • 10""^ reported by the NA31 collaboration at CERN indicated direct 
CP violation, the result of the E731 collaboration at Fermilab 0, (7.4±5.9) ■ 10~^, was still 
compatible with a vanishing value. The new measurement of the KTeV collaboration 0, 

Re{e'/€) = (28.0 ± 4.1) ■10-\ (1) 

is in agreement with the CERN experiment NA31 and rules out the superweak models. 
Additional information will be provided in the near future by the NA48 collaboration 
and by the KLOE experiment at DA$NE. In view of the new experimental result, whose 
statistical uncertainty will be further reduced in the future, it is particularly interesting to 
investigate whether the quoted range and the weighted average can be accommodated in 
the standard model. 

Direct CP violation measures the relative phases of the decay amplitudes for 

ir° -^ 7r°7r° and K^ -^ Tr+Tr". 

The two pions in these decays can be in two isospin states, / = (A/ = 1/2) and 1 = 2 
(A/ = 3/2). The two amplitudes acquire phases through final state strong interactions 
and also through the couplings of weak interactions. We can use Watson's theorem |^ to 
write them as 

{n7i,I\nw\K'^) = Aie'^', (2) 

{nn,I\nw\K'^) = A*je'^', (3) 

with 6i being a phase of strong origin which is extracted from tt — vr scattering. The 
remaining amplitude Aj contains a phase of weak origin. Throughout the paper we use 
the following isospin decomposition: 

A{K' ^ vr+TT-) = yf ^e^'° + ^^26^'^ , (4) 



A{K+ -^ 7r+7r°) = J^A^e''^ . (6) 



The parameter of direct CP violation can be written as 

e' ^ { Im^2 ImAo 



£ V2\e\ V ReAa ReA 



(7) 



with uj = ReA2/ReAo = 1/22.2. In Eq. (0) we used the fact that, to a large degree of 
accuracy, the strong interaction phases of e' and e cancel in the ratio (see e.g. Ref. 0). 
In order to obtain the numerical value of e' /e it is now necessary to calculate the two 
amplitudes (ImAo and Im742) including their weak phases. 

Using the operator product expansion, the K ^ titt amplitudes are obtained from the 
effective low-energy hamiltonian for \AS\ = 1 transitions [^ |], ^, 

din) = Zi{n) + ryiifi) , T = -Xt/Xu , Xq = V*, V^^ . (9) 

The arbitrary renormalization scale /i separates short- and long-distance contributions to 
the decay amplitudes. The Wilson coefficient functions Ci{^) contain all the information on 
heavy- mass scales. The terms with the ZiS contribute to the real parts of the amplitudes Aq 
and A2. The ?/j's, on the other hand, contribute to the imaginary parts and are relevant for 
CP violating processes. The coefficient functions can be calculated for a scale /i > 1 GeV 
using perturbative renormalization group techniques. They were computed in an extensive 
next-to-leading logarithm analysis by two groups ||T^, |TT[]. The Wilson coefficients depend 
on the CKM elements; the ?/j's are multiplied by At which introduces CP violation in the 
amplitudes. Finally, the calculation of the decays depends on the hadronic matrix elements 
of the local four-quark operators 

{Qi[^J))I = (tttt, /| g,(/i) |ir°) , (10) 

which constitute the non-perturbative part of the calculation. This is the main subject 
of this paper. The hadronic matrix elements will be calculated using the l/N^. expansion 
within the framework of the effective chiral lagrangian for pseudoscalar mesons [0, |1^, ^ . 
In a previous article |l^ we already reported the results of 0{p^/Nc) for the operators Qe 



and Qs- In this article we investigate one-loop corrections for all matrix elements relevant 

for e' /e. 



The local four-quark operators Qi{^) can be written, after Fierz reordering, in terms 
of color singlet quark bilinear s: 



Qs = 4 J^ SLl^dL qn^iiqR , 
g 

3 

Qt = 4 ^ -eg SLl^di qRlt^QR ■ 



Q4 = 4 ^ SL7''gL gL7^rfL : 



Qg 



-8 2^ -eg sxgij gijdL 



(11) 
(12) 

(13) 
(14) 



where the sum goes over the light flavors (g = m, d, s) and 



gii,L = 2(1^75)?: 



(2/3, -1/3, -1/3) 



(15) 



Qs - Q& arise from QCD penguin diagrams involving a virtual W and a c or t quark, with 
gluons connecting the virtual heavy quark to light quarks. They transform as (8^,, 1r) 
under SU{?>)l x SU{?>)r and contribute, in the isospin limit, only to A/ = 1/2 transitions. 



Q^ and Qs are electroweak penguin operators |T^, [T^. 

The imaginary parts of the amplitudes occurring in Eq. (|^ are those produced by the 

weak interaction. Thus we obtain the amplitudes 

r^ I 

ImAi = — ^ ImAi | ^2/i(/i) (Qi)/ • 



(16) 



Since the phase originating from the strong interactions is already extracted in Eq. (|^), 
absolute values for the J2i Vi {Qi)i should be taken. We shall return to this point later on. 
Collecting all terms together we arrive at the general expression 



e' 
e 



Gf 



UJ 



2 \e\ReAn 



ImAi 



n 



-n. 



112 

UJ 



with 



no = I X^ Viif^) {Qi)o 

i 

n2 = 1 5^ 2/^(/^) (Q^)2 



l-Q. 



■q+ri' j ) 



(17) 

(18) 
(19) 



where Q 



v+v' 



0.25 ±0.10 takes into account the effect of the isospin breaking in the quark 



masses (m„ 7^ nid) [|16], |T^ |T8[. We have written Eq. ([T7| ) as a product of factors in order to 



emphasize the importance and uncertainty associated with each of them. The first factor 
contains known parameters and takes the numerical value Gp 1^/(2 |£:| ReAo) = 346 GeV~^. 
The remaining terms are discussed in the following sections. Especially important to this 
analysis are the operators Qe and Qs which dominate the J = and 1 = 2 contributions in 
Eqs. (|T8|) and (pJ]), respectively. The terms i/i {Qi)2 are enhanced by the factor l/u, and a 
crucial issue is whether the enhancement is strong enough to produce an almost complete 
cancellation with the i/i {Qi)o terms, leading to an approximately vanishing e'/e even in 
the presence of direct CP violation, or whether the cancellation of the two terms is only 
moderate and a large value of e'/e can be obtained within the standard model. 

The paper is organized as follows. In Section ^ we briefly recall the numerical values of 
the CKM elements relevant to this analysis. In Section ^ we review the general framework 
of the effective low-energy calculation and discuss the matching of short- and long-distance 
contributions to the decay amplitudes. The next two sections contain our results, which 
we present in two steps. As the chiral theory is an expansion in momenta, we keep the 
first two terms in the expansion and calculate one-loop corrections to each term of the 
expansion separately. Loop corrections to the lowest terms, in the momentum expansion, 
are presented in Section ^ they are of OijP /N^ for the density operators and of OijP' /Nc) 
for the current operators. In Section ^ we extend the one-loop corrections to the next 
order in momentum by calculating corrections to the density operator Qq of 0{p'^/Nc). 
Numerical results for e'/e are included in both of these sections. Finally, our conclusions 
are contained in Section ||. 

2 The CKM Elements 

The second factor in Eq. (|17D originates from the CKM matrix. In the Wolfenstein para- 



metrization IS 



lm\ = Im(V;:v;,) = A^X'r] = KslK^I'^^, (20) 

since A = Vus and A = |Vcfe|/A^. Numerical values for the matrix elements are taken from 
the particle data group [Q and from Ref. [pl[] : 



\Vus\ = 0.2196 ±0.0023, (21) 

\Vcb\ = 0.040 ±0.002, and (22) 

\Vub\ = (3.56 ±0.56) ■ 10-^ (23) 

The last parameter we need is the phase 1] which is obtained from an analysis of the 
unitarity triangle whose overall scale is given by the value of \Vcb\- Such t] versus p plots 



are now standard [^, 23, 24, 25, 26| and are obtained, primarily, from |Kife/V"cb| which 
produces a circular ring and from a hyperbola defined from the theoretical formula for e. 
The position of the hyperbola depends on mt, iKbl, and I3k- The intersection of the two 
regions, together with constraints from the observed B^ — B^ mixing parameterized by 
AMd and the lower bound on B^ — B^ mixing, defines the physical ranges for t] and p. 



A very recent analysis can be found in Ref. [0. The remaining theoretical uncertainties 



in this analysis are the values of the non-perturbative parameters B^ in e, Fb^V Bd in 
(AM)d, and ^ = Fb.^/WsI FB^^/Bdm {AM)d/{AM)s. Bk has been calculated by various 
methods which, unfortunately, give a large range of values. Two recent calculations are 
found in Refs. 



Taking Bk = 0.80 ± 0.15, Fb^\/ Bd = 200 ± 40MeV 0, and 



28 



^ = 1.14 ± 0.08 [^ ^], the authors of Ref. |2^ obtain the following range for ImAj 



1.04 ■ 10"' < ImAi < 1.63 ■ 10" 



(24) 



where the experimentally measured values and the theoretical input parameters are scanned 
independently of each other, within the ranges given above. In Section ^ we shall use this 
range in the numerical analysis of e' /e. 

3 General Framework 



The method we use is the 1/Nc expansion introduced in Refs. |T^, 0. In this approach, 
we expand the hadronic matrix elements in powers of external momenta, p, and the ratio 
l/Nf.. In an earlier article WM we investigated one-loop corrections to lowest order in the 
chiral expansion for the operators Qq and Qs- The calculation of the one- loop corrections 
for current-current operators was done in Ref. |^, where predictions for the A J = 1/2 
rule were reported. 

To calculate the hadronic matrix elements we start from the effective chiral lagrangian 
for pseudoscalar mesons which involves an expansion in momenta where terms up to 0{p^) 
are included I^Sj. Keeping only (non-radiative) terms of 0(p'^) which are leading in Nc, 
for the lagrangian we obtain: 

f2 

eff = ^[{D.U^D^^U) + ^{\nU^ -\nU)^ + r{MU^ + UM^)) 



L 






(25) 



with DfJJ 



d,U 



ir^U + iUl^, {A) denoting the trace of A and M. = diag(m„, m^, m^). 
/^ and r^ are left- and right-handed gauge fields, respectively, / and r are parameters 
related to the pion decay constant F^ and to the quark condensate, with r = —2{qq)/ p. 
The complex matrix ?7 is a non-linear representation of the pseudoscalar meson nonet: 





[/ = exp -n , n = TT^Xa , 


where, in terms of the physical states. 




( 7," + -j^^arj + ^^hi y27r+ 


n = 


\/2t^~ -7r° + -^arj + 




i 72^- 727^0 



(AaAfe) = 2(5, 



ab ■ 



(26) 



hrj' 



V2K+ 



-73^^ + 



(27) 



3«^ / 



with 



cos^ — v2sin6' , 



V26 = sin^ + v^cos^. 



f28) 



The conventions and definitions we use are the same as those in Refs. ||I^ ^. In particular, 
we introduce the singlet tjq in the same way and with the same value for the Ua{^) symmetry 
breaking parameter, a = mi + mi, — 2m^ ^ 0.72 GeV , corresponding to the r] — 7]' mixing 



angle 9 = —19° |^J]. The bosonic representations of the quark currents and densities 
are defined in terms of (functional) derivatives of the chiral action and the lagrangian, 
respectively: 



ss 



qiLTQjL 



-i^-^{U^d^U)^. 



J« 



QiRljL 






-e// 



5 Mi 



■(—U^ + L^d^JJ^d^UU^ + 2rLsU^MU^ + rH2M^ 



(29) 
, (30) 

and the right-handed currents and densities are obtained by parity transformation. Eqs. 
(PPP and (^) allow us to express the four-fermion operators in terms of the pseudoscalar 
meson fields. The low-energy couplings Li, L2, and L3 do not occur in the mesonic densities 
in Eq. (|30|) . Furthermore, at tree level they do not contribute to the matrix elements of the 
current-current operators and have been omitted in Eq. (pUj). It is now straightforward to 
calculate the tree level (leading- A*";,) matrix elements from the mesonic form of the 4-quark 
operators. 

For the 1/Nc corrections to the matrix elements {Qi)i we calculated chiral loops as 
described in Refs. O, E9i. The factorizable contributions, on the one hand, refer to the 



strong sector of the theory and give corrections whose scale dependence is absorbed in 
the renormahzation of the effective chiral lagrangian. This property is obvious in the case 
of the (conserved) currents and was demonstrated exphcitly in the case of the bosonized 
densities [|14[ ^. Consequently, the factorizable loop corrections can be computed within 
dimensional regularization. The non-factorizable corrections, on the other hand, are UV 
divergent and must be matched to the short-distance part. They are regularized by a finite 



cutoff which is identified with the short-distance renormahzation scale [13, M, E^, M, R6l. 



The definition of the momenta in the loop diagrams, which are not momentum translation 



invariant, was discussed in detail in Refs. [|1J, ^. A consistent matching is obtained 
by considering the two currents or densities to be connected to each other through the 
exchange of a color singlet boson and by assigning the same momentum to it in the long- 



and short- distance regions |Tj, ^7\, ^, ^ ^, ^ . 

For the short- distance coefficient functions, we use both the leading logarithmic (LO) 
and the next-to-leading logarithmic (NLO) values.^ The published values for the Wilson 



coefficients are tabulated for scales equal to or larger than 1 GeV |T0|, |TT], |4^ . In Appendix ^ 



we present tables for the coefficient functions at scales 0.6 < fi < 1.0 GeV calculated with 
the same analytic formulas and communicated to us by M. Jamin [^. The NLO values are 
scheme dependent and are calculated within naive dimensional regularization (NDR) and 
in the 't Hooft-Veltman scheme (HV), respectively. The absence of any reference to the 
renormahzation scheme dependence in the effective low-energy calculation, at this stage. 



prevents a complete matching at the next-to-leading order p2[- Nevertheless, a comparison 
of the numerical results obtained from the LO and NLO coefficients is useful in order to 
estimate the uncertainties associated with it and to test the validity of perturbation theory. 

4 Analysis of e'/s 

In the twofold expansion in powers of external momenta and in 1/A^c we must investigate, at 
next-to-leading order, the tree level contributions from the 0{p'^) and the 0{p'^) lagrangian, 
and the one-loop contribution from the 0{p'^) lagrangian, that is to say, the 1/A'"c corrections 
at lowest order in the chiral expansion. In this section we combine our results and report 
values for e' /e up to these orders. 



^We treat the Wilson coefficients as leading order in 1/Nc since the large logarithms arising from the long 
renormahzation group evolution from (?Tit, Mw) to /i ~ 0(1 GeV) compensate for the 1/Nc suppression. 



4.1 Long-Distance Contributions 



The hadronic matrix elements for all the operators are calculated following the method 
described in the previous section. As mentioned above, we consider the bilinear quark op- 
erators to be connected to each other through the exchange of a color singlet boson, whose 
momentum is chosen to be the variable of integration. This is our matching procedure 
described in Ref. |jl^. For current- current operators, the tree level contributions from the 
0{p'^) and 0{p^) lagrangian and the one-loop contribution from the 0{p^) lagrangian are 
(9(p^), 0{p'^), and OijP' /Nc), respectively. For density-density operators they are OijP), 
OijP'), and 0{p^/Nc), respectively. The numerical results for the matrix elements to these 
orders are given in Tabs. |l| and |^ as a function of the cutoff scale. These values were 



obtained in Refs. |14, E9[ using the following values for the various parameters [20 



IT 


= (m^o+m^+)/2 = 


= 137.3 MeV, 


F^ -- 


= 92.4 MeV 


K 


= {rriKO +mK+)/'2 = 


= 495.7 MeV, 


Fk -- 


= 113 MeV 


V 


= 547.5 MeV , 




9 -- 


= -19°, 


V' 


= 957.8 MeV. 









Note that the matrix elements generally contain a non-vanishing imaginary part (cutoff 
independent at the one-loop level) which comes from the on-shell (tt — vr) rescattering. 

It is customary to parameterize the hadronic matrix elements in terms of the bag factors 
Bl and Bi \ which quantify the deviations from the values obtained in the vacuum 
saturation approximation |]4^ : 



B) 



(1/2) 



B. 



(3/2) 



Re(Qj)o 
Re(Q^)2 



I G {1,...,8}, 
I e {1,2,7,8}. 



(31) 
(32) 



The VSA expressions for the matrix elements are taken from Eqs. (XIX.ll)- (XIX.28) 
of Ref. ^^ , and the corresponding numerical values can be found in Refs. [0, ^ . We list 
the bag parameters in Tabs. ^ and |[0 One might note that the values of the B factors 
contain the real parts of the matrix elements and not their absolute values. For the 



amplitudes appearing in Eqs. (|T^ - (|1^) we need both real and imaginary parts for the 
matrix elements. We calculated the imaginary parts in the 1/A'"c expansion and included 

^The definition of the bag parameters in Eqs. ( |3l| ) and ( p^ ) refers to the complete sum of the factorizable 
and non-factorizable terms in the hadronic matrix elements. Therefore we are free of any of the infrared 
problems discussed in Ref. pa] , which occur for Qg with other definitions of i?g . 



A, 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Im((5i)o 


{Qi)o 


-33.2 


-40.2 


-48.2 


-57.3 


-67.4 


-5.55i 


{Q2)0 


58.8 


68.8 


79.9 


92.4 


106 


11. li 


{Q3)0 


0.05 


0.03 


-0.02 


-0.12 


-0.26 





(Q4)0 


92.1 


109 


128 


150 


173 


16.6i 


(^5)0 


-0.05 


-0.03 


0.02 


0.12 


0.26 





(Q6)0 


-38.6 


-33.7 


-29.4 


-25.5 


-21.9 





{Q7)0 


40.1 


46.6 


54.1 


62.6 


72.2 


8.32z 


{Q8)0 


119 


119 


119 


118 


117 


36. 7 i 



Table 1: Real and imaginary parts (last column) for the hadronic matrix elements of 
Qi,...,5,7 (in units of 10^ ■ MeV^) and Qe,8 [ in units of R"^ ■ MeV, with R = 2m^/(m, + m^) ] . 
The values are for the / = amplitudes in the isospin limit (m„ = nid) and for various 
values of the cutoff Ac. 



Ae 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


Im(g,)2 


(Ql)2 


2.51 


-2.26 


-7.77 


-14.0 


-21.1 


-3.45i 


(^2)2 


2.51 


-2.26 


-7.77 


-14.0 


-21.1 


-3.45i 


(^7)2 


-10.7 


-6.27 


-1.15 


4.67 


11.2 


5.18i 


(Q8)2 


35.3 


31.2 


27.2 


23.2 


18.8 


-n.5i 



Table 2: Same results as in Tab. |T] for the 1 = 2 amplitudes. 



their values in Tabs. [T| and |^. They are produced by the imaginary parts of the one-loop 
diagrams, as required by unitarity. In order to study the sensitivity of the results on the 
imaginary part we calculated the matrix elements by two methods (for a discussion of this 
point see Ref. |2^). In the first method, we obtain the absolute values by correcting the 
real parts using the phenomenological phases. This procedure has also been followed in 
Refs. 



, 46 



In the second method we assume zero phases and use the real parts of 
the matrix elements. The second method, to a large degree of accuracy, corresponds to 
adopting the phases obtained in the 1/A'"c expansion. 



Analytical formulas for all matrix elements are given in Refs. p4| , p9|l . Among them 



10 



Ac 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


5(1/2) 


8.24 


9.98 


12.0 


14.2 


16.6 


4^/^) 


2.91 


3.41 


3.96 


4.57 


5.23 


5(1/^) 


0.004 


0.002 


-0.002 


-0.010 


-0.021 


5f/'^ 


2.54 


3.00 


3.53 


4.13 


4.75 


5(1/2) 


0.0009 


0.0005 


-0.0003 


-0.0014 


-0.0020 


5(V^) 


1.10 


0.96 


0.84 


0.72 


0.62 


5}^/^) 


0.16 


0.18 


0.21 


0.23 


0.26 


S(V2) 


1.21 


1.21 


1.21 


1.20 


1.19 



Table 3: Bag parameters for the / = amplitudes, shown for various values of the cutoff. 
i?5 7 8 depend on i? ~ 2m\/ms and are calculated for a running ras{^ = Ac) at the leading 
logarithmic order (Aqcd = 325 Me V) with ms{l GeV) = 175 MeV. 



Ae 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


5(3/2) 


0.11 


-0.10 


-0.34 


-0.61 


-0.92 


i?f^) 


0.11 


-0.10 


-0.34 


-0.61 


-0.92 


Bf"^ 


-0.10 


-0.06 


-0.01 


0.04 


0.09 


r(3/2) 


0.64 


0.56 


0.49 


0.42 


0.34 



Table 4: Same results as in Tab. H for the 1 = 2 amplitudes. 



four are particularly interesting and important, and we repeat them here: 



(Qi)c 



(Q2)C 



-T/' (' 



nip^ — 771, 



I) 



. 4Lg 2 1 

1 + -^^mt + 



F2 - (4^)2ir^2 



X (6A2 
2 



c V ^'nnj^ + Qm^)\ogK^ + 



v^ 



F^ {ml - ml) 



1 4L^ 2 1 

1 + -^ml + 



F2 



(4vr)2F,2 



15 



11 2 _ 15 
^ 4 



X ( ^A^ + ( ^mi - ^ml ] log A2 + 



(33) 



(34) 



11 



(Q6)0 

{Qs)2 



4 



AVS / 2m 

F^ \ms + rrid 



[m 



2ml 



2^2 \ms + md 
1 



K 



1 + 



m^ 



jr 



16(47r 



;^ logAe + 



(35) 



im, 



K 



F3 



a -2a]-^{3Li- %u 



F3 



(47r)2F2 



m, 



K 



wi„ 



g-ajlogAc 



(36) 



where the elhpses denote finite terms, which for brevity are not written exphcitly here, but 
have been included in the numerical analysis (in particular, they provide the mass terms 
which make the logarithms dimensionless) . The constants Lg and Lg are renormalized 
couplings and defined through the relations 0] 



'K 



F^ 



1 + 



4L^ 
F3 



{m\ - ml) 



(37) 



and 



m 



K 



m + m. 



mt 



1- 



[m, 



K ''■kJ / f r o f^\ 
— 1^^5 - ZLf^) 



(3J 



Their values are Lg 



2m 

2.07 ■ 10^^ and Lg = 1.09 • 10^^. Lg has a small dependence on the 
ratio m^s/rrid and we shall use the (central) value (of) m^s/rrid = 24.4 ± 1.5 ||47[| . 

In Eqs. (^) - (|36|) we have summed the factorizable contributions [first two terms on 
the r.h.s. of Eqs. (P^) and (|5^), first term on the r.h.s. of Eq. (^), and first three terms on 
the r.h.s. of Eq. (^) ] and the non-factorizable contributions. Factorizable terms originate 
from tree level diagrams or from one-loop corrections to a single current or density whose 
scale dependence is absorbed in the renormalization of the effective chiral lagrangian (i.e., 
in Ftt, Fk, L^, Lg, and in the renormalization of the masses and wave functions). Finite 
terms from the factorizable loop diagrams for {Qq)o and (^8)2 are not absorbed completely 
and must be included in the numerical analysis [|T^. We note that the couplings Li, L2, 
and L3 do not contribute to the matrix elements of Qq and Qs to O(p^) and to 0{p'^) for the 
current-current operators. The non-factorizable contributions, on the other hand, are UV 
divergent and must be matched to the short-distance part. As we already discussed above, 
they are regularized by a finite cutoff, Ac, which is identified with the renormalization scale 
/i of QCD. 

We discuss next Eqs. (^3]) - ( P^ ) which have several interesting properties []T1|, ^|. First, 
the VSA values for {Qi)o and {Q2)o are far too small to account for the large AJ = 1/2 
enhancement observed in the CP conserving amplitudes. Using the large- iVc limit [ -BJ = 
3.05, B2 = 1.22] improves the agreement between theory and the experimental result, 
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but it still provides a gross underestimate. However, the non-factorizable l/N^. corrections 
in Eqs. (^) and (P^ ) contain quadratically divergent terms which are not suppressed with 
respect to the tree level contribution, since they bring in a factor of A = A^/(47rF7r)^ and 
have large prefactors which, in some cases, can be as large as six in Eq. (^). Quadratic 
terms in (Qi)o and {Q2)o produce a large enhancement (see Tab. H) which brings the AJ = 
1/2 amplitude in agreement with the observed value p9[. Corrections beyond the chiral 
limit {niq = 0) in Eqs. (|3^) - (|3^) are suppressed by a factor oi 6 = m\- ^ / {An Ft,Y ~ 20 % 
and were found to be numerically small. 

The case of {Qq)^ and {Qs)2 is different from that of (Qi,2)o- The leading-A^c values 
are very close to the corresponding VSA values. Moreover, the non-factorizable loop cor- 
rections in Eqs. (^) and (PB|), which are of 0{p^/Nc), are found to be only logarithmically 
divergent [l^. Consequently, in the case of {Q8)2 they are suppressed by a factor of 6 



compared to the leading C(p°) term and are expected to be of the order of 20 % to 50 % 
depending on the prefactors. We note that Eq. p6| ) is a full leading plus next-to-leading 
order analysis of the Qg niatrix element. The case oi Bq is more complicated since the 
0{p^) term vanishes for Qq. Nevertheless, the non-factorizable loop corrections to this 
term remain and have to be matched to the short- distance part of the amplitudes. These 
0{p^/Nc) non-factorizable corrections must be considered at the same level, in the twofold 
expansion, as the 0{p'^) tree contribution. Consequently, a value of Bq around one 
[which corresponds to the 0{p'^) term alone] is not a priori expected. However, numer- 
ically it turns out that the 0{jP /N^ contribution is only moderate (see Tab. ^. This 
property can be understood from the {U'^)dq{U)qs structure of the Qq operator which van- 
ishes to 0{p^) implying that the factorizable and non-factorizable 0{p^ /N^ contributions 
cancel to a large extent [1^. The fact that the factorizable and non-factorizable terms to 



this order have infrared divergences which must cancel in the sum of both contributions 



gives another qualitative hint for a value oi B\ remaining around one [Q. This explains 



why for Qg to 0{p^/Nc) we do not observe a A/ = 1/2 enhancement similar to the one for 

(-]/2) fS/2~) 

Qi and Q2 in the CP conserving amplitude. The leading- A'"c values for Bq and Bg 
are therefore more efficiently protected from possible large l/N^ corrections of the 0{p'^) 

lagrangian than 5| 2 . The effect of the 0{p'^ /Nc) term is however important for B^ as 

(3/2) 
for B\ because it gives rise, in general, to a noticeable dependence on the cutoff scale 

I^ , which is relevant for the matching with the short- distance part (see below). We note 



that i?g shows a scale dependence which is very similar to the one of Sg (see Tabs. ^ 
and ^ leading to a stable ratio i?g /-Bg over a large range of scales around the value 

(A Icy) /Q /o~j 

Bq V-Bg — 1.77 ± 0.05 where the error refers to the variation of Ac between 600 MeV 
and 1 GeV. The 0{p^ /N^) corrections consequently make the cancellation of Qg and Qg in 
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e' je less effective. 

Finally, B\ and B2 were found to be too small to account for the measured value 
of the A/ = 3/2 amplitude for small values of the cutoff and even become negative for large 
values of Ac (see Tab. ^. Due to an almost complete cancellation of the two numerically 
leading terms [^ B\ and B\ are expected to be sensitive to corrections from higher 
order terms and/or higher resonances.^ For this reason, even though the effect of 5} 2 
is small, in the analysis of e' je we will not use the values listed in the table but we will 
extract the parameters B\ and B2 from the experimental value of Re742. This point 
is further discussed in the next section. 



4.2 Numerical Results 

Collecting together the values of the matrix elements in Tabs. 1 and 2 and the values of 
the Wilson coefficients in Appendix ^ we can give now the numerical results for e' je. As 
we already mentioned above, we use the real parts of the matrix elements and consider two 
cases. In the first case, we use the real part of our calculation and the phenomenologically 
determined values for the final state interaction phases, (5o = (34.2 ± 2.2)° and 82 = 
(-6.9 ± 0.2)° gH, in order to arrive from Eqs. (|1|) and (|T9|) to 

Ho = ^- Vi/.(/i)Re(g,)o(l-fi,+,0, (39) 

COSOo "^ 

n2 = ^5^2/.(^)Re(Q,)2. (40) 

COS02 "^ 

The factor l/cos6j enhances the A/ = 1/2 term in Eq. (^) by about 25 % with respect to 
the A/ = 3/2 one and consequently makes the cancellation between the Qq and Qg opera- 
tors even less effective. It allows us to estimate the effect of multiple (vr — vr) rescattering 
on the imaginary part. In the second case, we assume zero phases and use the equations: 

no = 5^1/.(yt^)Re(g,)o(l-fi,+,'), (41) 

i 

Ha = 5^i/.(/i)Re(g,)2. (42) 



•^As explained in Ref. ||2^, the sensitivity of i?} 2 to corrections from higher order terms is expected 
to be smaller. Therefore, the fact that the l/Nc expansion, at this stage, does not reproduce the A/ — 3/2 
amplitude does not imply that the A/ = 1/2 amplitude cannot be calculated to a sufficient degree of 
accuracy. This point was also illustrated in Ref. ||2^ where it was shown that higher order corrections 
investigated with a Nambu Jona-Lasinio model arc much larger for the A/ = 3/2 channel than for the 
A/ = 1/2 one. 
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The comparison of the two cases provides, in part, an estimate for higher order effects. The 
latter case gives numerical results very close to those we would get if we used the imaginary 
parts from Tabs. |l] and |^. As we already mentioned, we extract the values of i?| and 
i?2 from the experimental value for ReA2. This procedure has also been followed in the 
phenomenological approach of the Munich group (last reference of |I0| ) . Then 

R/n\r^ R/n\r^ v^ cos (52 ReA^ 8.42 ■ 10^ MeV=^ 

Re Qi 2 /i = Re Q2 2 /i = ^ .. .. / , w X = / , w X , 43 

with ReA2 = 1.5 ■ 10"^ MeV. The values of zi and Z2, for 600 MeV < /i < 1 GeV, are listed 
in Appendix A of Ref. |2^. All other B factors are taken from Tabs, ^and ^ In particular, 
for B[ and B2 we use the values listed in Tab. ^. These numbers were obtained in 



Ref. p9[ where it was shown that they saturate the observed value of KbAq and are in good 



agreement with the phenomenological result of Ref. fl^.^ We note that Qi and Q2 do not 



give a direct contribution to e'/e since yi and ?/2 are zero. Rather, {Qi)o,2 and ((52)0,2 are 
used to sum up the contributions from the operators (^4, Qg, and Qio which are redundant 
below the charm threshold (see below). 

The elimination of the scale dependence of QCD in the numerical result is an impor- 
tant criterion and we discuss it in some detail. The Wilson coefficients in the effective 
hamiltonian in Eq. (^ depend on the renormalization scale /i. This should be matched 
with the scale dependence of the chiral operators and their respective matrix elements. 
The bosonization of the density-density operators introduces masses which are also scale 
dependent. In particular, ((56)o and ((53)2 are proportional to i?^ = [2m|./(ms + mrf)] ~ 
Am'^/m'j. which brings in a /i dependence through the quark masses already for the tree 
level (factorizable) contributions. This is different from the matrix elements of the current- 
current operators which are /i independent in the large- A''^ limit. In the products of y^ and 
ys with the corresponding matrix elements, the fi dependence from the running quark 
mass is exactly cancelled by the diagonal evolution of the Wilson coefficients taken in the 
large- A'c limit |]T^ ^. This property is preserved at the two- loop level [0. Furthermore, 
the yU dependence beyond the m^ evolution, i.e., the /i dependence oi Bq and Bg 
was shown in QCD to be only very weak for values above 1 GeV [|l^. This requires that 
the (non-factorizable) l/Nf. corrections to the matrix elements of the Qq and Qg operators 
(which produce the scale dependence of the B factors) should not show a large dependence 



"^Even though not fully consistent from a theoretical point of view, the values for S} 2 i'^ Tab. ^ can 
be used together with the experimental value for Re^o in the prefactor of Eq. ([l7[), since the numbers in 
the table produce a 

e'/e is rather small. 



the table produce a value for Re^o close to the experimental one [^ . In addition, the effect of 5} 2 ^^ 
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on the cutoff scale. The fact that the 0{p^/Nc) terms in Eqs. (|35|) and ( p6D have only 
a logarithmic cutoff dependence is for this reason welcome. Finally, the decrease of both 
B factors with Ac = /i in Tabs. ^ and |^ which is due to these logarithms is qualitatively 



consistent with their /i dependence found for /i > 1 GeV in Ref. ||T^, i.e., it has the correct 
slope. As shown below the residual scale dependence of -Bg and Bg even if moderate 
is still too large to allow an exactly scale independent result for e' /e. 

Throughout the numerical analysis of direct CP violation we take Aqcq = A^ = 
325 ± 80MeV corresponding roughly to as{Mz) = 0.118 ± 0.005. For fi^+r;' we adopt 
the range given after Eq. (plQl). The status of the strange quark mass has been reviewed 



recently in Refs. |22|, |2^ , and we use the range 

m, (IGeV) = 150 ± 25 MeV, (44) 

which is in the ball park of the values obtained in the quenched lattice calculations (see 
Ref. |5D| and references therein; for a very recent analysis see Ref. [^]) and from QCD sum 



rules [Q. We note that the QCD sum rule results are generally higher than the lattice 
values. Lower bounds on the strange mass have been derived in Ref. [p^ . 

In Fig. we depict e' /e as a function of the matching scale (/i = A^), calculated from 
Eqs. ( p9D and (^OD with LO Wilson coefficients for the central value of ImAt and for various 
values of m^, flri+ri', and Aqcd according to their ranges defined above. For low values of 
the matching scale we find a rather moderate enhancement of the VSA result which is due 
to the weaker cancellation between the Qq and Qg operators. However, one might note 
that very large values for e'/e in the range of the recent Fermilab measurement 0] are not 
reached with the B factors listed in Tabs. ^ and ^ together with Eq. (^3D, if central values 
are used for the parameters. Indeed, adopting central values for m^, fi^+^z, Aqcd, and ImA* 
and varying A^ between 600 MeV and 900 MeV we obtain as 'central range' for the CP 
ratio: 

8.4 ■ 10^^ < e'/e (LO-central) < 12.9 ■ 10"^ . (45) 

This is also illustrated in Fig. |^ where we show the various contributions to e'/e for central 
values of the parameters at a scale of Ac = 700 MeV. For this value of the cutoff B^ ' 
is very close to unity whereas i?g is significantly suppressed which leads to a value for 
e'/e of 11.5 ■ 10~^. Smaller numbers are obtained for larger values of the cutoff. Another 
noticeable contribution, beside that of Qs, which reduces the value of e'/e is the / = 
component of Qi and Q2. As we already mentioned above, this contribution comes from 
the Q4, Qg, and Qio operators which are redundant below the charm threshold and satisfy, 
to LO and at NLO in the HV scheme, the relations in Eq. (^6]). In the NDR scheme the 
relations receive small 0{as) and 0{a) corrections ||10|| . 
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Figure 1: e'/e using LO Wilson coefficients and the experimental phases, plotted for various 
values of ms(l GeV) and 1^,,+,,' as a function of the matching scale Ac = /i. We use ImA^ = 



1.33- 10 '^. The sohd (dashed, dot-dashed) lines correspond to A 



QCD 



325 (245, 405) MeV. 



In Fig. ^we show how the various terms depend on the choice of the matching scale. In 
particular, we observe that the behaviour of e'/e is almost identical to the one of y^ (Q6)o- 
This is due to the fact that the ratio Bq /B^ is approximately stable over the whole 
range of the cutoff Ac and, consequently, i/q {Qg)q and (— ) ys {Q%)2 fall off roughly in the 
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Figure 2: The J = and 1 = 2 contributions to e'/e of the dominant operators using LO 
Wilson coefficients and the experimental phases. We also use the central values for Aqcd, 
ms(l GeV), O^n+r]', and ImAj at a scale of Ac = 700 MeV. The contributions of the operators 
Qs and Qs are negligible and are not included in the figure. 

same way. We note that the ratio yQ{fi)/ml{fi) increases by about 12 % if the scale /x = Ac 
is varied between 600 MeV and 1 GeV, whereas B^ (Ac) decreases by 44%. A similar 
statement applies to Qg. The decrease of i?g and i?g is therefore qualitatively con- 
sistent with the (non- diagonal) evolution of i/q and ys computed in the leading logarithmic 
approximation, and it leads to a fairly moderate overall scale dependence. This property 
is due to the fact that the 0{p^ /N^ terms in Eqs. (|35|) and (|36|) have only a logarithmic 
cutoff dependence which, nevertheless, still goes beyond the /i dependence of the short- 
distance part. In this situation it would be tempting to adopt the large- A^c values for B^ 

(3/2) 

and -Bg which are scale independent and coincide, to a very good approximation, with 



their VSA values B, 



(1/2) 



B. 



(3/2) 



1. However, the results show that 1/Nc corrections 
are important, and to recover the VSA values would require an a priori unexpected can- 
cellation of the 0{p^/Nc) corrections with higher order terms or contributions from higher 
resonances. Therefore, the VSA might underestimate the true range of uncertainty in the 
analysis of e'/e. 
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Figure 3: Various contributions to e'/e using LO Wilson coefficients and the experimental 
phases, plotted as functions of the matching scale Ac. We use central values for m^, fir]+ri', 
Aqcd, and ImAi. 



The dependences of the result on ms, ^-q+r)', and Aqcd are given in Fig. |1|. Among 
them the m^ dependence is the most important one. The dependence on ImAj, to a large 



degree of accuracy |27], is multiplicative and can be obtained in straightforward way. If 
we take into account the residual dependence on the matching scale by varying /i = A^ 
between 600 and 900 MeV and scan independently the theoretical input parameters and 
the experimentally measured numbers, we obtain the following range for e'/e calculated 
with LO Wilson coefficients: 

3.1 -10'^ < e'/e (LO) < 31.6 -10-^ (46) 

The quoted range results from a variation of m^, fi.^+.^/, and Aqcd in Eqs. (|39|) and (0) as 
depicted in Fig. |1| and also allows for a variation of ImA^ according to the range defined 
in Eq. (|2|). 

We investigate next the dependence on the NLO Wilson coefficients. The NLO values 
are scheme dependent and are calculated within naive dimensional regularization (NDR) 
and in the 't Hooft-Veltman scheme (HV), respectively. As already mentioned, the absence 
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Figure 4: Same as in Fig. [1|, now using NLO Wilson coefficients in the NDR scheme. The 
sohd (dashed, dot-dashed) hues correspond to Aqcd = A^ = 325 (245, 405) MeV. 



of any reference to the renormahzation scheme dependence in the low-energy calculation 



prevents a complete matching at the next-to-leading order p2[. Nevertheless, a comparison 
of the numerical results obtained from the LO and NLO coefficients is useful in order to 
estimate the corresponding uncertainties and to test the validity of perturbation theory. 
In the NDR scheme, introducing the NLO coefficients does not noticeably affect our 
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Figure 5: Same as in Fig. ^, now using NLO Wilson coefficients in the NDR scheme. 





Case 1 


Case 2 


LO 


8.4 < e'/e < 12.9 


6.3 < e'/e < 9.5 


NDR 


8.0 < e'/e < 11.0 


5.9 < e'/e < 8.4 


HV 


5.8 < e'/e < 6.6 


4.2 < e'/e < 4.7 


LO + NDR + HV 


5.8 < e'/e < 12.9 


4.2 < e'/e < 9.5 



Table 5: Central ranges for e'/e (in units of 10~^) at LO and NLO (NDR and HV). 
The numbers are obtained for central values of rus, fi,y+,y', ImA^, and Aqqd by varying A^, 
between 600 and 900 MeV. 'Case 1' and 'Case 2' correspond to the use of Eqs. (H)- (|0D 
and (|4l|) - (^21) , respectively. 



numerical results (see Fig. H). For A^ = Aqqd ^ 325 MeV we find slightly lower values for 
e'/e and a somewhat larger difference between the results obtained for Aqcq = 325 MeV and 
405 MeV, respectively. Generally, the difference between LO and NLO is more pronounced 
for very low values of the matching scale, but it is still moderate except for Aqcd = 405 MeV. 
For Aqcd = 325 MeV (245 MeV) the effect of the NLO coefficients is rather small, and values 
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Figure 6: Same as in Fig. Rl, now using NLO Wilson coefficients in the HV scheme. 





Case 1 


Case 2 


LO 


3.1 < e'/e < 31.6 


2.4 < e'/e < 23.2 


NDR 


2.7 < e'/e < 26.4 


2.1 < e'/e < 20.2 


HV 


1.9 < e'/e < 16.5 


1.5 < e'/e < 11.9 


LO + NDR + HV 


1.9 < e'/e < 31.6 


1.5 < e'/e < 23.2 



Table 6: Same results as in Tab. ^ but for the complete scanning of the parameters (A^, 
TUs, ^n+ri' and Aqcq, and ImAt) as explained in the text. 

for the matching scale as low as 600 -700 MeV appear to be acceptable. We also notice 
a slightly smaller scale dependence, that is to say, the NLO Wilson coefficients further 
improve the stability of the calculation. This property becomes obvious if we investigate 
the various contributions to e'/e (compare Figs. § and §). In particular, at NLO in the 
NDR scheme we observe a smaller variation of i/q {Qq)o and ys {Q8)2 in the range of Ac 
between 600 MeV and 1 GeV. Nevertheless, the numerical effect of the NLO coefficients is 
rather moderate, and the 'central' and scanned ranges quoted in Tabs. || and ^ are close 
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Figure 7: Same as in Fig. 0, now using NLO Wilson coefficients in the HV scheme. The 

QCD — ^4-MS 



sohd (dashed, dot-dashed) hues correspond to A, 



Ag^ = 325 (245, 405) MeV. 



to the LO results given in Eqs. (^) and 

In the HV scheme, the effect of the NLO coefficients is more pronounced. Both y^ {Qq)o 
and ys {Qs)2 are rather stable over a large range of the matching scale leading to an 
approximately stable result for e' /e between 700 MeV and 1 GeV. This is shown in Fig. ^ 
for the central values of ms(l GeV), ^r]+r]', IniAi, and Aqcd- On the other hand, for Aqcd = 
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405 MeV and A^ < 700 MeV we observe a noticeable slope indicating the breakdown of the 
perturbative expansion of QCD (see Fig. |^. However, for moderate values of the matching 
scale the numerical values for the ratio depend weakly on the choice of Aqcd (see Fig. |^), 
which makes the result rather stable. Generally, at NLO in the HV scheme we obtain 
smaller values for e' je (see Tabs. |] and |^). 

We note that at NLO the maximum value for the ratio is found for moderate values 
of Ac around 700 - 800 MeV, whereas the upper bound in Eq. (^61) refers to low values of 
the matching scale (~ 600 MeV). Finally, one might note that the numerical values of the 
Wilson coefficients in the HV scheme communicated to us by M. Jamin [Q correspond to 
the treatment of the two-loop anomalous dimensions used in Ref. [|^ which differs from 
the one used in Ref. (ITTI. For this reason the NLO corrections to the Wilson coefficients 



in the HV scheme presented in Appendix ^ are generally smaller than the ones found 
in Ref. lTl| (for a discussion of this point see also Ref. ||2^ ). 

So far in the numerical analysis we have used Eqs. ( |39D - (^0|) together with the phe- 
nomenological values for the phases [^. Replacing them by Eqs. (^1]) - (|4^) leads to lower 
values for e' je (see Case 2 in Tabs. ^ and P). The numerical results are very close to those 
we would get if we used the imaginary parts obtained at the one-loop level in the 1/A^c 
approach. Both the central values and the upper bounds in the scanned ranges for e' je are 
lower due to a smaller contribution from the A/ = 1/2 terms. However, the modifications 
do not change substantially our picture of e' je. As mentioned above, the comparison of 
the two cases provides, in part, an estimate for higher order effects. 

In conclusion, the fact that we use rather low values for the matching scale makes 
some of the Wilson coefficients rather sensitive to NLO corrections. In particular, y^ and 
?/8 depend noticeably on the choice of the 75 scheme in dimensional regularization. For 
example, for ^ = 700 MeV and Aqcq = 325 MeV the values of y^ and ys at LO and in 
the NDR and HV schemes differ approximately by 20-30%. Since the non-perturbative 
calculation of the matrix elements is insensitive to this dependence, the corresponding 
uncertainty must be included in the final result for e' /e. Collecting together the LO and 
NLO values in Tab. || from Eqs. (^ - (gD|) and (|T|) - (g2|) we get the following range: 

4.2 ■ 10-^ < e'/e (central) < 12.9 ■ 10"^ , (47) 

which, for central values of rris, ^2,,+,,', Aqcq, and ImAi, takes into account the theoretical 
errors inherent to the method (dependence on the scheme and matching scale). Further- 
more, it includes the expected errors due to the neglect of higher order corrections to the 
imaginary part. Similarly collecting the values in Tab. ^ we obtain the following range 
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from the complete scanning of the parameters: 

1.5-10"^ < ^je < 31.6 ■10"\ (48) 

which also takes into account the uncertainties in the values for m^, H^+^z, Aqcd, and 
IniAi-l The upper bound from our calculation in Eq. (|48D is rather close to the central 
value of the new Fermilab measurement 0] and requires a conspiracy of the parameters 
within their ranges of uncertainties given above. 

The present world average for the ratio including earlier measurements is Re(£'/e) = 
(21.8±3.0)-10~^ 1^. Our result indicates that the experimental data can be accommodated 
in the standard model. A major uncertainty in the theoretical estimate of e' je is due to 
the choice of ttIs, which enters the calculation through the matrix elements of the operators 
Qe and Qg [see Eqs. (0) and (||)]. In Eq. (||) we have taken m,(l GeV) = 150 ± 25 MeV 



which is in the range of the values obtained in quenched lattice calculations and from the 
QCD sum rules. Adopting even lower values for m^ would allow us to relax the upper bound 
quoted above. However, recently the ALEPH collaboration analyzed the measured mass 



spectra of the strange r decay modes and reported a value oirrisiX GeV) = 234 1*1 ^g MeV [^ 



It will be interesting to see whether this large (central) value for vtis will remain when the 
error is reduced. Very recently, the value 7715(1 GeV) = (188 ± 22) MeV was obtained 
using a r-like decay sum rule for the meson |Q, which is consistent with the range 
used in this paper. The determination of ImAt will be further improved by precision tests 
of the unitarity triangle |^2[ removing to a large extent the corresponding uncertainty in 
the analysis of direct CP violation, fi^+^z which measures the contribution to e' je from 
the isospin breaking in the quark masses was estimated in Ref. |jl6| in the large- A'"c limit, 
and it will be a challenge to investigate, in future studies, the 1/A^c corrections to this 
parameter. Finally, the calculation of the hadronic matrix elements even though largely 
improved by including 1/A'"c corrections may still be plagued by noticeable uncertainties. 
Our analysis so far included terms of 0[jP\ 0{p^ /N^, and 0{p'^) for the matrix elements 
of the density-density operators and terms of 0{p'^), 0{p'^/Nc), and 0(p^) for the matrix 
elements of the current-current operators. In the following section we shall investigate the 
effect of higher order corrections. In particular we will consider the terms of 0{p'^ /N^ for 
the matrix elements of Qq. 



^A comparison with the results of other calculations performed within a specific scheme and treatment 
of the final state interaction phases should be done using the numbers in Tab. 0. 
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5 Higher Order Corrections 

In the previous section we have shown that the calculation of the hadronic matrix elements 
in the l/N^ expansion leads to a well defined range of values for e' /e which can account, to a 
large extent, for the weighted average of the experimental measurements 0, ^, ^. However, 
the central values obtained are lower than the values of the new measurement @]. The 
upper bound from our calculation requires, within the standard model, specific values of 
various parameters. In particular, lower values of the strange quark mass are favoured. In 
our analysis so far we varied the theoretical input parameters independent of each other 
and considered the experimental results within one standard deviation. This conservative 
attitude may to some extent exaggerate the differences [Q. In the present section we 
investigate the higher order corrections and consider in particular the question whether 
these corrections are able to substantially enhance the prediction for e' /e, so that a large 
value for the ratio could be explained even for central values of the input parameters. 

In the twofold expansion, the higher order corrections to the matrix elements of Qg 
and Qs are of orders: C(p^), 0{p^/N^), and OijP' /Nc). In this section we will consider 
the 0{p'^/Nc) contribution which brings in, for the first time, quadratic corrections on the 
cutoff. From general counting arguments we show that these corrections are expected to be 
large for Qq, which is a peculiar operator. Qq is consequently not protected from possible 
large corrections beyond the large- A''^ limit, and we cannot exclude the possibility that 
the contribution of Qq brings e' /e close to the experimental value for central values of the 
parameters. Calculating the 0{p^ /Nc) correction for Q^ in the chiral limit we explicitly 
find that it is indeed large and positive. 

Before investigating the 0{p^ /N^) corrections we briefiy estimate part of the higher 
order corrections replacing the '1/F^ expansion' by a 'I/Fr- expansion'. As already dis- 
cussed in Ref. [|T4[, we could have used the ratio 1// in place of 1/ F^^ in the next-to- leading 



order terms of Eqs. (^)- (|5§). This choice would be consistent at the level of first order 
corrections in the twofold series expansion, as the difference concerns higher order effects. 
However, the scale dependence of / (which is mainly quadratic) is absorbed through the 
factorizable loops to the matrix elements at the next order in the parameter expansion and 
does not occur in the matching with the short- distance contribution [0. Consequently, 



it is more appropriate to choose the physical decay constant in the expressions under con- 
sideration. In this situation, we can use, instead of F^r, the kaon decay constant F^ which 
gives an indirect estimate of higher order corrections. In Tab. ^ we show the effect of this 
modification, to 0{p^/Nc), on the values of -Bg and -Bg . We notice that the numbers 
are generally larger for the ''l/F^ expansion'. In particular, for A^ < 900 MeV the -Bg 
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Ac 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


r(1/2) 


1.10 


0.96 


0.84 


0.72 


0.62 




(1.30) 


(1.19) 


(1.09) 


(0.99) 


(0.91) 


^(3/2) 


0.64 


0.56 


0.49 


0.42 


0.34 




(0.72) 


(0.66) 


(0.59) 


(0.53) 


(0.46) 



>(l/2) 



?(3/2) 



Table 7: Bag parameters Bq and B^ shown for various values of Ac. The numbers in 
the parentheses are obtained by replacing F.,, by Fk in the next-to-leading order expression. 



factor is enhanced compared to the VSA value. This change, in spite of the somewhat 
smaller reduction of -Bg , leads to a moderate enhancement of e'/e which further im- 
proves the agreement with the observed value. Numerically, collecting together the various 
terms we get 5.3 ■ lO"'^ < e'/e (central) < 15.8 ■ 10~^. Scanning independently the input 
parameters we obtain 1.8 ■ 10~^ < e'/e < 38.3 ■ 10""^ in place of Eq. (^8]). Adopting zero 
phases reduces the upper bound to 28.3 ■ 10~^. 

In Fig. H we depict e'/e as a function of the matching scale (/x = A^) , calculated with 
LO and NLO (NDR and HV) Wilson coefficients and for the central values of ImAj and 
Q^ri+r)' and various values of mg. The variation of fi.^+,,' does not change the qualitative 
behaviour; it only shifts the curves upward or downward for smaller or larger values of 
fin+Tj', respectively. The curves in Fig. | result from replacing l/F,, by l/F^ in all next- 
to-leading order expressions relevant to the complete set of matrix elements. Beside the 
enhancement of the numerical result we also observe a somewhat smaller dependence on 
the matching scale. Finally, even though we obtain somewhat larger values for e'/e the 
effect is still rather moderate and does not affect the statement we made above that lower 
values of the strange quark mass are favoured. 

In the above we have argued that estimating the effect of higher order corrections 
to the matrix elements, by replacing the 'l/F^r expansion' by a '1/Fk expansion', does 
not drastically modify the results we obtained in the previous section. In particular, 
this statement refers to terms of 0{p^/N^) which are corrections on top of the 0{p^/Nc) 
contribution and correspond to the same pseudoscalar representation of the four-quark 
operator. In the following we will not study the 1/N^ corrections, which correspond to a 
two-loop diagram in the chiral theory. The same approximation was made in the chiral 



quark model ^g] . However, estimating the typical effect of higher orders by modifying the 
known corrections (l/i^jr —>' ^/ Fk) does not account for possible contributions from new 
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Figure 8: e' /e using LO and NLO (NDR and HV) Wilson coefficients and the experimental 
phases, plotted for various values of ing as a function of the matching scale Ac = /i, now 
with 1/-FV -^ ^/ Fk in the next-to-leading order terms for the matrix elements. We use 
the central values for ImAf and fi,j+r;'. The solid (dashed, dot-dashed) lines correspond to 
A^^^ = 325 (245, 405) MeV. 

terms which are absent at the level of the ffist order corrections. In particular, higher order 
terms in the p^ expansion (tree level) cannot be calculated because the low-energy couplings 
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in the 0{p^) lagrangian are very uncertain or even unknown. Nevertheless, these terms 
are independent of the (non-factorizable) matching scale and are chirally suppressed with 
respect to the leading term (~ m|^^/A^, with A^ ~ 1 GeV the scale of chiral symmetry 
breaking). The convergence of the tree level series was verified for the current- current 
operators |T^, ^ and also for Qg 10' where a complete leading plus next-to-leading order 



calculation now exists and tree contributions appear to decrease monotonically. For the 
operator Qq the leading term is 0{p'^) and by analogy we expect the higher order tree 
terms to be smaller. The terms of 0{p'^/Nc) on the other hand are not expected to be 
small for Qq. We remind the readers that for the CP conserving amplitude it is mainly 
the (quadratic) 0{p^ /N^) corrections which bring to {Qi^2)o a large enhancement relative 
to the (leading-A^c) C'(p^) values. As the leading-A^c value for Qq is also C(p^) we cannot 
a priori exclude that the value of {Q6)o is largely affected by OijP' jNc) corrections too. As 
already discussed in Section ^]l|, quadratic 0{p'^/Nc) corrections are proportional to the 
factor A = AI/IAttFt^Y relative to the O(p^) tree level contribution. Different is the case of 
the operator Qg since its leading- A'"c value is 0{p^) at lowest order in the chiral expansion. 
Quadratic terms for Qq are consequently chirally suppressed with respect to the leading- A^c 
value. More precisely the suppression factor is ~ (m|-_^/A^) ■ A. In contrast to {Qe)o, it is 
very unlikely that the 0{p'^/Nc) corrections for {Qs)2 could be larger than the 0{p^/Nc) 
contributions investigated in the previous section. 

We calculate next the 0{p'^/Nc) quadratic corrections to the matrix elements of the 
operator Qq. The pseudoscalar representation of Qq can be read off from Eq. (|30|) : 



Qe = -2fr'Yl 



\f{U^)d,iU),, + {U^)d,{L,Ud,U^d^U + 2rLsUM^U 



+rH2M)^^ + {L.U^d^Ud^U^ + 2rL^U^MU^ + rH2M^)^^{U),s 



+ 0{p^). (49) 



In the following we will calculate the 0{p'^/Nc) evolution of the operator Qq in the chiral 
limit. It is then straightforward to compute the hadronic matrix element {Qq)o. To calcu- 



late the evolution of Qq we use the background field method as described in Ref. |39] and 



also in Refs. |T^, 0. This operatorial method is very convenient to calculate corrections 
in the chiral limit. To this end we decompose the matrix U in the classical field U and the 
quantum fluctuation C,, 

U = exp{iV2^/f)U, e = r^, (50) 



with U satisfying the equation of motion 



a 



Nc 



Ud'^U^ -d'^UU^ + rUM^ -rMU^ = —{\nU-\nU^)-l, f/ = exp(i7r"A„//) . (51) 
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The O(p^) lagrangian thus reads 
C = C 



^(^/.r^^^a) +\{[d,^, ^wm^) - -^{euM^ + u^eM) 



aee+o{e). (52) 



The corresponding expansion of the meson density in Eq. 
yields 



around the classical field 



{D 



L)ij 



QiRQjL 

r 



n 



2 2 

+ i2V2 jLs [ W{^, MW} ],, +2^Ls[ W{e, MW} 



jji 



ij« 



Oie) ■ (53) 



Using Eq. (|53D the evolution of Qe can be obtained in a straightforward way. Integrating 
over the quantum fiuctuation by calculating the non-factorizable diagrams of Fig. |^.a we 
get the following result: 

3 A? 



Q,{Al) = -AFyilid'^U^d.U) 



ds 



1 + 



(54) 



This result has already been presented in Ref. |^ 
effect of the quadratic term in Eq. 



2{A7tF^)\ 
Before investigating the numerical 



a few comments are necessary: 

In Eq. (|5^) we present only the diagonal evolution, i.e., the term proportional to 
the operator {d^Wdf^U)ds which gives the only non-vanishing contribution to the 
K -^ TTTT amplitudes. This property is analogous to the tree level. One might note 
in particular that the Lg contribution vanishes since it does not produce a term 
proportional to this operator. The H2 contribution vanishes from the beginning as 
it does not appear in Eq. (|5c 



To 0{p^/Nc) we showed explicitly that the factorizable contributions provide the 
corrections needed to obtain the physical values of the low-energy couplings [^ 



Except for finite corrections, the values of the couplings can be obtained in the large- 
Nc limit, i.e., by imposing tree level relations in order to set up the renormalized 
(factorizable) matrix elements [compare Eqs. (0) and (p8|)]. To 0{p'^/Nc) in Eq. ( p^ 
we use again the renormalized coupling Lg, defined in the large- A^c liniit|, since the 



^Note that our constants L[ should not be confused with the rcnormahzation scale dependent coefficients 
ir inRefs. M and M. 
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Figure 9: Non-factorizable loop diagrams for the evolution of Qq to 0{p'^ /N^). The crossed 
circles denote the two bosonized densities, the black circles the strong interaction vertices 
from the kinetic term in the lagrangian. Similar diagrams with the 0{p^) mass or Ua{^) 
terms in Eq. (|52| ) are logarithmically divergent or finite, and we do not present them here. 

scale dependence of the bare coefficient L5 will be absorbed in factorizable loop 
corrections to the matrix elements and does not need to be matched to any short- 
distance contribution WM.. 



Beside the diagrams in Fig. |^.a a priori the diagram in Fig. |^.b with a strong vertex 
proportional to L5 could also contribute. However, it is easy to see that since the 
L5 term in the lagrangian contains a quark mass matrix M., this diagram produces 
an operator similar to the one resulting from the Lg term at tree level for Q^. This 
operator does not contribute to the K ^ tttt amplitudes. 

The diagram of Fig. ^.b with the strong vertex coming from the Li, L2, L3, and Lg 
couplings also turns out to vanish. 



There are no 0{p'^/Nc) tree level contributions, i.e., from couplings L^ [33] which are 
subleading in Nc. 



Numerically, we observe a large positive correction from the quadratic term of 0{p'^ /Nc) 
in Eq. (^l)). The slope of this correction is qualitatively consistent and welcome since it 
compensates for the logarithmic decrease at 0{p^/Nc). Adding the 0{p'^/Nc) term to the 
full OijP') and 0{p^/Nc) result in Eq. (|35|) we obtain the following matrix element for Qq: 



(Q6)0 



4^3 



R\ml 



m^ 



i + ? 



A2 

c 



2(47r)2F^V 16(47r 



log A' + 



(55) 



The corresponding values oi Bq (obtained by adding the quadratic corrections to the 
values in Tab. |^) are listed in Tab. ^. The B^ factor is found to be rather stable around 
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Ac 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


r(1/2) 


1.50 


1.51 


1.55 


1.62 


1.73 



?(V2) 



Table 8: Bq ' , now including, in the chiral limit, terms of 0{p /N^ 



the value B, 



(1/2) 



1.6 ± 0.1. The quadratic term of 0{p'^/Nc) is of the same magnitude 



as the 0{p'^) tree term. For the (Qi)o and ((52)0 there was also a large enhancement to 
the 0{p'^) tree contribution introduced by the 0{p'^/Nc) term (see Tabs. ^ and ^. Qq is a. 
A/ = 1/2 operator, and the enhancement of (Q6)o suggests that at the level of the l/Nc 
corrections the dynamics of the A J = 1/2 rule also applies to Qq. One might however 
note that the long-distance evolution of the operator Qq including both the OijP /Nc) and 
0{p'^/Nc) terms is very different from the one of Qi or Q2. The former is approximately 
constant over a wide range of the cutoff scale due to the smaller coefficient of the quadratic 
term and a large cancellation of the scale dependences between the quadratic and the 
logarithmic term, whereas the later [which does not receive any 0{p^ /N^) contribution] has 
a large positive slope. One should also remark that we observe a noticeable suppression of 
{Qs)2 similar to the one needed for (Qi,2)2 in order to bring the CP conserving A/ = 3/2 
amplitude down to the experimental value.Q 

Using the values for Sg in Tab. ^ together with the bag factors of the remaining 
operators presented in the previous section we calculated again the ratio e' /e for central 



values of r/io, Vt 



s 1 ^ '"q+Tj' 5 



ImA/, and A 



QCD- 



The results for the three sets of Wilson coefficients 



LO, NDR, and HV and for Ac between 600 and 900 MeV are given in Tab. p. The numbers 
are obtained from Eqs. (|39| ) - (^0|) and (|4l|) - (|42|) , respectively. With LO Wilson coefficients. 



the enhancement oi Bq leads to larger values for e' /e, and the predictions are now more 
stable and closer to the data. The results in the NDR scheme are rather close to the LO 
ones although more sensitive to the value of Aqcd- In the HV scheme, the effect of the 
NLO coefficients is more pronounced. The results are significantly smaller for low values 
of the matching scale and less stable. This is illustrated in Fig. [1^ where we show e'/e for 
various values of mg as a function of the matching scale. 

Performing a scanning of the input parameters as explained above, we arrive at the 
values in Tab. |10[ Comparing these results with the values in Tabs. ^ and ^ we see a clear 
enhancement originating from the quadratic term in Eq. (|55D . The large ranges reported 



^The mechanism for the suppression of {Qb)2 however differs from the one for (Qi,2)2, since in the 
former case it occurs through logarithms and in the latter case mainly through quadratic terms Jo, p9[. 
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Case 1 


Case 2 


LO 


19.5 < e'je < 24.7 


14.8 < e'/e < 19.4 


NDR 


16.1 < e' je < 23.4 


12.5 < e'/e < 18.3 


HV 


9.3 < e'/e < 19.3 


7.0 < e'/e < 14.9 


LO + NDR + HV 


9.3 < e'/e < 24.7 


7.0 < e'/e < 19.4 



Table 9: Same results as in Tab. ^ but now including quadratic terms of 0{p'^/Nc) for Qq 
as explained in the text. 





Case 1 


Case 2 


LO 


8.0 < e'/e < 62.1 


6.1 < e'/e < 48.5 


NDR 


6.8 < e'/e < 63.9 


5.2 < e'/e < 49.8 


HV 


2.8 < e'/e < 49.8 


2.2 < e'/e < 38.5 


LO + NDR + HV 


2.8 < e'/e < 63.9 


2.2 < e'/e < 49.8 



Table 10: Same results as in Tab. |^, but now including quadratic terms of 0{p^/Nc) for 
Qq as explained in the text. 

in Tab. [T^ can be traced back to the large ranges of the input parameters. The parameters, 
to a large extent, act multiplicatively, and the larger range for e'/e is due to the fact that 
the central value(s) for the ratio are enhanced roughly by a factor of two compared to the 
results we presented in the previous section. More accurate information on the parameters, 
from theory and experiment, will restrict the values for the CP ratio. 

The contributions from current-current operators to e'/e are rather small, and correc- 
tions from higher order terms (beyond the ones investigated in Section ^) and from higher 
resonances will not be able to modify their B factors in a way that they change the ratio 
e'/e considerably. For the reason explained earlier in this article, higher order corrections 
for the operator Qs are not expected to bring a large change for the ratio. The question of 
whether a large value of e'/e can be accommodated in the standard model without specific 
values of various parameters reduces essentially to the value of -Bg . Adding the 0{p'^/Nc) 
quadratic terms produces a substantial increase for the value of the matrix element {Q6)o, 
and a large value for B^ in the range of 1.6 cannot be excluded. This property leads 
to a more natural explanation for a large value of e'/e. Our result can be modified by 
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Figure 10: e'/e using LO and NLO (NDR and HV) Wilson coefficients and the experimental 
phases, plotted for various values of rris as a function of the matching scale Ac = /i, including 
for Qq the term of 0{p^ /N^) in the chiral limit. We use again the values ImAj = 1.33 ■ 10~^ 
and Vtr^+r)' = 0.25. The solid (dashed, dot-dashed) lines are for Aqcd = 325 (245, 405) MeV. 



corrections of 0{p'^/Nc) beyond the chiral limit, from logarithms and finite terms, but they 
are not expected to remove the large enhancement observed in Eq. (p5|). Nevertheless it 
would be very interesting to verify this statement through an explicit calculation. 
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In view of the large corrections one might question the convergence of the 1/A^c ex- 
pansion. However, we remind the readers that the quadratic term of 0{p'^/Nc) in Qq we 
consider in this section represents a new class of terms absent to 0{p'^) and 0{p^ /N^). It is 
reasonable to assume that the terms of 0{p'^/Nc) (~ A) carry a large fraction of the entire 
(non-factorizable) contribution, since quadratic corrections in the cutoff from higher order 
terms are chirally suppressed (i.e., they are ~ A ■ (5, A ■ 6"^, . . .). 

We point out that the quadratic terms obtained at the level of the pseudoscalar mesons 
are physical and must be included in the numerical analysis. Our result is compatible 
with the fact that, in a complete theory of mesons, the quadratic dependence on the cutoff 
should be absent. Indeed one expects that incorporating higher resonances allows one to 
select higher values for the cutoff and does not remove the effect of the quadratic terms, but 
turns them smoothly into logarithms. Therefore, within a limited range of the cutoff, the 
quadratic terms provide an approximate representation of the effect of higher resonances. 
This behaviour has been observed in the calculation of the vr^ — 7r° mass difference after 



including the vector and axial vector mesons [57, 58|. In this particular case the quadratic 



terms turn into finite terms. It has also been observed partly for the B^ parameter after 



including the vector mesons [^. The two examples in Refs. |5^, ^ show that effects 
of higher resonances can modify noticeably the results and reduce the dependence on 
the cutoff but are not strong enough to reverse the effects of the pseudoscalar mesons. 
Nevertheless, it would be very interesting to include higher resonances in the calculation 
of the K ^ nir decays, in order to study explicitly their effect. 

We note that the A/ = 1/2 enhancement we observe for Qe to 0{p'^/Nc) is not able 
to render the contribution of Qq in KeA^ dominant. The numbers for Bq in Tab. ^ 
are consistent with the observed value for the A/ =1/2 amplitude, which is dominated 
by Qi and Q2', at LO and at NLO in the HV scheme, the enhancement of the Bq 
factor changes the result of Ref. |^ by less than 5% for ms(lMeV) = 150 MeV and 
245 MeV < Aqcd < 405 MeV. In the NDR scheme, the effect can amount to approximately 
11 % of the amplitude for Aqcd = 405 MeV. Therefore we do not see a correlation between 
the large values of ReAg and e' /e, since the two quantities are dominated by different 
operators.0 In particular, in one of the first estimates of e'/e 0] it was suggested (following 



Ref. [§]) that the A J = 1/2 amplitude is dominated by the operator Qq, which would lead 
to a large value of e'/e. Such a mechanism, in the range for the matching scale we consider, 
would require an enhancement of {Qe)o several times larger than the one obtained in the 
present analysis. The same remark applies to Ref. []61 |. 



^For a discussion of this point see also Ref. |21 
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Among the previous calculations, loop corrections to the operators Qq and Qs in the 
l/Nf. expansion were also considered in Refs. P3| , |3B[. This study used a different matching 
condition, and the parametrization of the 0{p^) lagrangian was not general. The authors 
obtained a large reduction of (Q8)2 and an enhancement of {Qq)q, predicting large values 
for e' /e. 

It is interesting to compare our values for e' /e with the results found with other methods. 
Lattice calculations obtained Bf''^\2 GeV) = 1.0 ± 0.2 and Bf''^\2 GeV) = 1.0 ± 0.2 and 
predicted a small value for e' /e = (4.6±3.0) ■ 10~^ with Gaussian errors for the experimental 
input (see Ref. ||62| and references therein). More recent values reported for B^ are 
5f/^^(2GeV) = 0.81(3)(3) H, 0.77(4)(4) H, and 1.03(3) ||. 5^^^^^ was estimated 
in Ref. [§§: 5^^/^^(2GeV) = 0.76(3)(5). However, as stressed in Ref. ^, the systematic 



uncertainties in this calculation are not completely under control. This statement has been 
confirmed by a recent analysis ^^ which obtained negative values for Bq and favours 



either negative or slightly positive values for e'/e. Although the scales used in lattice 
calculations and the phenomenological approaches are different, the various results for the 
B factors can be compared, for values of the scale around 1 GeV or above, since B^ 

(3/2) 

and Bg were shown in QCD to depend only very weakly on the renormalization scale 



for values above 1 GeV [jTO[- Small values for e'/e consistent with zero were also quoted 



in Ref. 68 



The chiral quark model pUf yields a range for i?g in the HV scheme which is above 



6 

the VSA value, Bg (0.8 GeV) = 1.6 ± 0.3, and predicts a small reduction of the B^ 
factor, 5^/^^(0.8 GeV) = 0.92 ± 0.02. The quoted range for the CP ratio is 7 ■ 10"^ < 
e'/e < 31 ■ 10^^. Since the treatment of the renormalization scale in Ref. |^ is different 



from the one used in this article we do not see a clear link which could easily explain why 
both approaches give approximately the same result for Bq . 

Very recently, an extensive study of e'/e in the standard model was presented in 
Ref. [^. The authors investigated the sensitivity of the CP ratio on the input param- 
eters and updated their numerical values. They treated B^ and i?g as parameters 

(1/2) (3/2) 

and adopted the values B^' = 1.0 ± 0.3 and -Bg = 0.8 ± 0.2 together with the con- 
straint B^^^^^ > Bf''^\ Numerically, they obtained 1.05 ■ 10"^ < e'/e < 28.8 ■ 10"^ and 
0.26 ■ 10"^ < e'/e < 22.0 ■ 10"^ in the NDR and HV schemes, respectively. The quoted 
results are consistent with the values we get for e'/e to OijP') and OijP /N^ for Qq and 
Qs- As shown in this section, an additional large contribution comes from the 0{p'^/Nc) 
term of the Qq operator. 

In summary, we have shown in this section that the quadratic terms of 0{p^ /Nc) are 
large for Qq. From general counting arguments we have good indications that among the 
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various next-to-leading order terms in the p^ and l/N^. expansions they are the dominant 
ones. They enhance i?g and bring e' /e much closer to the measured value for central 
values of the input parameters. We obtain a quadratic evolution for Qq which indicates 
that a A J = 1/2 enhancement mechanism is operative for Qq as for (5i,2- -Bg i^ expected 
to be affected much less by terms of 0{p^ /Nc) due to an extra p^ suppression factor relative 
to the leading 0{p^) tree term. 

One should recall that our analysis is performed in the chiral limit. Corrections be- 
yond the chiral limit, from logarithms and finite terms, are not expected to remove the 
large enhancement of i?g arising from the quadratic term in Eq. (pS]).^ Consequently, 
the results for e' /e we obtained in Section ^, by including terms of OijP') and 0{p^/Nc) 
for Qq and Qs, should be considered as a lower range, which is shifted to higher values 
by including also quadratic terms of OijP' jN^. The ideal case would be to calculate and 
include the full OijP' jNc) amplitudes, as well as the 0{j)'^) and 0{jP /N^) terms. It would 
also be interesting to investigate the effect of higher resonances (at least the vector mesons 
and presumably also the axial vector and scalar mesons). Each of the additional effects 
separately is not expected to counteract largely the enhancement found for B^ . Nev- 
ertheless, in the extreme (and unlikely) case where all these effects would come with the 
same sign a significant modification of the result cannot be excluded formally. In order 
to reduce the scheme dependence in the result for e' je appropriate subtractions would be 
necessary (see Refs. [^ WW)- 

In view of the noticeable uncertainties connected still with both the calculation of the 
matrix elements and the exact values of the various parameters taken from theory and 
experiment, it is difficult to decide whether the large value of e' /e observed recently is 
indicating new physics beyond the standard model [^ |69|, ^. In this situation it is 
interesting to investigate other kaon decays in order to perform precision tests of flavour 
dynamics and to search for new physics [|22|, [7T| . 

6 Conclusions 

In this article we have presented results of a new analysis for the CP parameter e' /e. 
Our interest in this topic concentrates on the improved calculation of loop corrections for 
the hadronic matrix elements. It is well known that the leading-iVc values for the matrix 
elements underestimate the A/ = 1/2 amplitude Aqvh K ^ tttt decays. It has been shown 

^Note that for the CP conserving A/ = 1/2 amphtude the chiral Umit gives a good approximation to 
the numerical result. 
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earlier that l/N^ contributions to the operators Qi and Q2 are very large, bringing the 
value for the amplitude Aq closer to the experimental value [0; an improved matching 
condition brings it even closer [^]. The same method introduced corrections to the matrix 



elements of the operators Qe and Qs [0, HI |36| and modified the predictions for the 
parameter e' /e. 

In view of this knowledge and the fact that three large experiments [0, |^, ^ were 
measuring the CP asymmetry, we decided to embark on an extensive study of the hadronic 
matrix elements including at the same time improvements of the input parameters which 
have taken place in the meantime. In particular, we incorporate an improved estimate 
of the multiplicative CKM factor [^] and use leading and next-to-leading order Wilson 
coefficients, which were communicated to us by M. Jamin [^]. 

In the first part of the paper (up to Section 4) we have presented our results to C(p^) 
and 0{jP /Nc) for the dominant operators Q^ and Q^ fl^ and have included them in an 
extensive analysis of the CP parameter. We have found that the matrix elements {Q%)q 
and {Qs)2, to this order, have only a logarithmic dependence on the cutoff. The corrections 
to these operators are smaller than those of {Qi)o and ((^2)0 which are quadratic in the 
cutoff 1^. They decrease {Q8)2 roughly to half its value in the VSA and modify {Qq)o 
to a lesser extent leading to a ratio Bq /B^ ~ 1.8. The net effect is to eliminate 
the almost complete cancellation between the two operators but the overall values of the 
matrix elements are reduced. The corresponding ranges for e' /e are given in Tabs. |^-^ 
and Figs. |l]-|^. Adopting central values for the input parameters (m^, ^2^+^/, A^, and 
ImAi) we obtain 4.2 ■ 10"'^ < e'/e (central) < 12.9 ■ 10"'*. The quoted range refers to the 
uncertainties associated with the calculation of the hadronic matrix elements and with the 
use of three sets of Wilson coefficients LO, NDR, and HV. Performing a complete scanning 
of the parameter space we obtain 

1.5-10"* < e'/e < 31.6- 10~* . 

The upper values for e'/e obtained in this part of the article are close to the experimental 
data 1^, ^. They are reached only for low values of rris and specific values for the other 
parameters. As stated in the article, this is the complete first-order calculation for Qq in 
the twofold expansion and provides a benchmark for additional corrections. 

A major part of the present article is the estimate of still higher order effects. In this 
direction we have studied, first, the changes introduced by the replacement of the coupling 
constant F-^ by Fk in the next-to-leading order expressions for the matrix elements, which 



gives an indirect estimate of higher order corrections |I^. We found that the predicted val- 
ues are increased. Numerical results for central values of ImAf and Qr]+r]' and various values 
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of rUs are shown in Fig. |, which indicate that the experimental data can be accommodated 
in the standard modeL A low value of m^ is also favoured. 

In a second step we studied the 0{p'^/Nc) corrections for Qq. Here the 0(p°) term 
vanishes; the 0{p^/Nc) correction was found to be moderate |l^. Thus a significant cor- 
rection appears, for the first time, through quadratic terms of 0{p'^/Nc), and the behaviour 
of (Q6)o is similar to the one found for the matrix elements {Qi)o and {Q2)o [iHI- In Sec- 
tion 5 we have shown that the value for {Qq)o is enhanced by the 0{p^ /N^ contribution 
in the chiral limit. This point we already emphasized in Ref. [Q. Numerically, we ob- 
tain values for B\ around 1.6 ± 0.1. Our calculation indicates that at the level of the 
1/N(. corrections a A/ = 1/2 enhancement is operative for Qq similar to the one of Qi 
and Q2 which dominate the CP conserving amplitude. The effect of adding the 0{p'^/Nc) 
quadratic terms is evident as a substantial increase in the value of e' /e, which brings the 
result rather close to the data for central values of the input parameters. Numerically, this 
is shown by the following range obtained by collecting together the results for the three 
sets of Wilson coefficients LO, NDR, and HV: 

7.0-10-^ < e'/e (central) < 24.7-10-^ 

(for details see Tab. |^). Performing a complete scanning of the parameter space for the 



various cases produces the ranges reported in Tab. ITU. 

As stated in the article, it is still desirable to calculate and compare the full amplitudes 
to 0{p^), 0{p'^/Nc), and 0{p^/N^). The incorporation of higher resonances would be very 
interesting since it would allow to select higher values for the matching scale. A more 
sophisticated treatment of the scheme dependence remains a challenge for future studies. 
However, it is encouraging that the approximations we made in this paper give results 
close to the experimental data. Clearly the possibility of an natural explanation, within 
the standard model, of the experimental value for e'/e cannot be excluded. 

To sum up, we have presented our results from an extensive study of the hadronic 
matrix elements to 0{p'^/Nc). We computed all matrix elements in the same theoretical 
framework, except for {Qi)2 = (^2)2 which were extracted from the data on the CP 
conserving decays. Our predictions for e'/e are close to the weighted experimental average 
for central values of the input parameters. 

Note added: After completion of this article the NA48 collaboration at CERN reported 
the value Re (e'/e) = (18.5 ± 7.3) ■ 10""^ |]72|. The new world average is Re (e'/e) = 



(21.2 ± 4.6) ■ 10 ^. The conclusions of the present article are in agreement with this new 
measurement. 
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A Numerical Values of the Wilson Coefficients 



In this appendix we list the numerical values of the LO and NLO (HV and NDR) Wilson 
coefficients for AS = 1 transitions used in Section 14.21 These values were communicated 



to us by M. Jamin [^. Following the lines of Ref. [T^ the coefficients i/i are given for a 10- 
dimensional operator basis {Qi, . . . , Qio}- Below the charm threshold the set of operators 
reduces to seven linearly independent operators [see Eqs. ( pA]) - (0)] with 



Q4 



-Qi + Q2 + Q3 



Q9 = IQi - IQs , Qio = IQi + Q2 - ^Qs . (56) 



At next-to-leading logarithmic order in (renormalization group improved) perturbation 
theory in the NDR scheme the relations in Eq. (|56D receive 0{as) and 0{a) corrections 
jlOl , ^ . In the present analysis we use the linear dependence at the level of the matrix 
elements {Qi)i, i.e., at the level of the pseudoscalar representation where modifications to 
the relations in Eq. (|56|) are absent. We note that the effect of the different treatment of 
the operator relations at next-to-leading logarithmic order, which is due to the fact that in 
the long-distance part there is no (perturbative) counting in a^, is numerically negligible. 
The following parameters are used for the calculation of the Wilson coefficients: 



M, 



w 



80.2 GeV, 



sin^^ 



w 



0.23, 



mt 



170 GeV, 



mblrrib 



4.4 GeV, 



a 



mArric 



1/129, 



1.3 GeV. 



/^ 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


1/2 


0.0 


0.0 


0.0 


0.0 


0.0 


ys 


0.0410 


0.038 


0.035 


0.034 


0.032 


1/4 


-0.056 


-0.056 


-0.055 


-0.055 


-0.055 


1/5 


0.009 


0.011 


0.012 


0.012 


0.013 


ye 


-0.133 


-0.116 


-0.106 


-0.098 


-0.092 


yi/a 


0.024 


0.025 


0.027 


0.028 


0.029 


ys/a 


0.217 


0.180 


0.155 


0.138 


0.125 


2/9 /a 


-1.749 


-1.657 


-1.595 


-1.550 


-1.515 


yio/a 


1.007 


0.887 


0.803 


0.740 


0.690 



Table 11: A5 = 1 LO Wilson coefficients for A 



QCD 



245 MeV. 
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/i 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


Vs 


0.052 


0.046 


0.043 


0.040 


0.038 


Vi 


-0.063 


-0.063 


-0.062 


-0.062 


-0.062 


y5 


0.008 


0.010 


0.012 


0.013 


0.013 


y& 


-0.187 


-0.154 


-0.135 


-0.122 


-0.113 


yi/a 


0.029 


0.031 


0.033 


0.034 


0.036 


1/8 /« 


0.324 


0.249 


0.206 


0.178 


0.158 


1/9/a 


-1.957 


-1.799 


-1.702 


-1.634 


-1.585 


1/10/tt 


1.280 


1.082 


0.956 


0.867 


0.800 



Table 12: A5 = 1 LO Wilson coefficients for A 



QCD 



325 MeV. 



/x 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


yz 


0.065 


0.057 


0.051 


0.048 


0.045 


1/4 


-0.069 


-0.070 


-0.069 


-0.069 


-0.069 


1/5 


0.005 


0.009 


0.011 


0.013 


0.014 


y& 


-0.285 


-0.209 


-0.173 


-0.152 


-0.137 


1/7/tt 


0.033 


0.035 


0.038 


0.039 


0.041 


ys/tt 


0.526 


0.356 


0.277 


0.230 


0.198 


1/9/a 


-2.295 


-1.995 


-1.836 


-1.736 


-1.666 


l/io/tt 


1.690 


1.334 


1.139 


1.011 


0.920 



Table 13: AS* = 1 LO Wilson coefficients for A 



QCD 



405 MeV. 
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/i 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


Vs 


0.033 


0.031 


0.029 


0.028 


0.027 


Vi 


-0.050 


-0.051 


-0.051 


-0.051 


-0.050 


y5 


-0.014 


-0.005 


-0.001 


0.002 


0.004 


y& 


-0.154 


-0.120 


-0.103 


-0.093 


-0.085 


yi/a 


-0.037 


-0.035 


-0.034 


-0.033 


-0.032 


1/8 /« 


0.234 


0.189 


0.163 


0.146 


0.134 


1/9/a 


-1.783 


-1.658 


-1.586 


-1.537 


-1.502 


1/10/tt 


0.971 


0.803 


0.700 


0.630 


0.577 



Table 14: AS = 1 NLO Wilson coefficients (NDR) for A 



QCD 



A. 



(4) 



245 MeV. 



/x 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


yz 


0.040 


0.036 


0.034 


0.032 


0.030 


1/4 


-0.055 


-0.054 


-0.053 


-0.053 


-0.053 


1/5 


0.017 


0.015 


0.014 


0.014 


0.014 


y& 


-0.125 


-0.103 


-0.090 


-0.082 


-0.077 


1/7/tt 


-0.033 


-0.033 


-0.033 


-0.032 


-0.032 


ys/tt 


0.269 


0.212 


0.181 


0.160 


0.146 


1/9/a 


-1.791 


-1.663 


-1.588 


-1.539 


-1.503 


l/io/tt 


0.990 


0.816 


0.711 


0.638 


0.585 



Table 15: AS = 1 NLO Wilson coefficients (HV) for A, 



QCD 



A. 



(4) 



245 MeV. 
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/i 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


Vs 


0.037 


0.038 


0.036 


0.034 


0.032 


Vi 


-0.051 


-0.056 


-0.057 


-0.058 


-0.058 


y5 


-0.067 


-0.024 


-0.011 


-0.004 


-0.001 


y& 


-0.334 


-0.199 


-0.150 


-0.126 


-0.111 


yi/a 


-0.052 


-0.037 


-0.034 


-0.032 


-0.031 


1/8 /« 


0.413 


0.289 


0.229 


0.195 


0.173 


1/9/a 


-2.160 


-1.864 


-1.718 


-1.633 


-1.576 


1/10/tt 


1.445 


1.079 


0.889 


0.771 


0.690 



Table 16: AS = 1 NLO Wilson coefficients (NDR) for A 



QCD 



A. 



(4) 



325 MeV. 



/x 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


yz 


0.052 


0.048 


0.043 


0.040 


0.037 


1/4 


-0.064 


-0.063 


-0.062 


-0.061 


-0.061 


1/5 


0.037 


0.021 


0.018 


0.017 


0.016 


y& 


-0.229 


-0.155 


-0.124 


-0.107 


-0.097 


1/7/tt 


-0.023 


-0.030 


-0.031 


-0.031 


-0.030 


ys/tt 


0.500 


0.329 


0.255 


0.214 


0.188 


1/9/a 


-2.188 


-1.875 


-1.724 


-1.636 


-1.577 


l/io/tt 


1.489 


1.102 


0.904 


0.783 


0.699 



Table 17: A5 = 1 NLO Wilson coefficients (HV) for A, 



QCD 



A. 



(4) 



325 MeV. 



44 



/i 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


Vs 


-0.020 


0.039 


0.043 


0.041 


0.039 


Vi 


-0.012 


-0.056 


-0.063 


-0.065 


-0.065 


y5 


-0.447 


-0.092 


-0.036 


-0.017 


-0.008 


y& 


-1.327 


-0.415 


-0.244 


-0.181 


-0.149 


yi/a 


-0.218 


-0.054 


-0.036 


-0.032 


-0.031 


1/8 /« 


0.788 


0.488 


0.342 


0.269 


0.227 


1/9/a 


-3.019 


-2.236 


-1.927 


-1.768 


-1.672 


1/10/tt 


2.422 


1.538 


1.162 


0.958 


0.829 



Table 18: AS = 1 NLO Wilson coefficients (NDR) for A 



QCD 



A. 



(4) 



405 MeV. 



/x 


0.6 GeV 


0.7 GeV 


0.8 GeV 


0.9 GeV 


1.0 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


0.0 


y2 


0.0 


0.0 


0.0 


0.0 


0.0 


yz 


0.023 


0.059 


0.055 


0.050 


0.046 


1/4 


-0.053 


-0.072 


-0.072 


-0.071 


-0.070 


1/5 


0.217 


0.048 


0.026 


0.021 


0.019 


y& 


-0.620 


-0.272 


-0.183 


-0.145 


-0.125 


1/7/tt 


0.051 


-0.019 


-0.026 


-0.028 


-0.028 


ys/tt 


1.206 


0.582 


0.385 


0.296 


0.246 


1/9/a 


-3.154 


-2.267 


-1.939 


-1.774 


-1.676 


l/io/tt 


2.589 


1.587 


1.187 


0.975 


0.842 



Table 19: AS = 1 NLO Wilson coefficients (HV) for A, 



QCD 



A. 



(4) 



405 MeV. 
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